Rare events and scaling properties in field-induced anomalous dynamics 



R. 



(N 

o 

(N 
O 

Q 



<Z3 



o3 



i 

o 
o 



> 

(N 
(N 



X 



Burioni, 1 G. Gradcnigo, 2 A. Sarracino, 2 A. Vczzani, 3 ' 4 and A. Vulpiani 5 

1 Dipartimento di Fisica e Scienza della Terra, Universita, degli Studi 
di Parma and INFN, viale G.P.Usberti 7/ A, 43126 Parma, Italy 
2 ISC-CNR and Dipartimento di Fisica, Universita Sapienza, p.le A. Moro 2, 00185 Roma, Italy 
3 Centra S3, CNR-Istituto di Nanoscienze, Via Campi 213A, 41125 Modena Italy 
'* Dipartimento di Fisica e Scienza della Terra, Universita degli Studi di Parma, viale G.P.Usberti 7/A, 43126 Parma, Italy 
5 Dipartimento di Fisica, Universita Sapienza and ISC-CNR, p.le A. Moro 2, 00185 Roma, Italy 

(Dated: December 7, 2012) 

We prove that, in particular regimes of a broad class of continuous time random walks, a small 
external field can turn diffusion from standard into anomalous. We illustrate our findings in the 
trap model, a prototype of subdiffusion in disordered and glassy materials, and in the Levy walk 
process, which describes superdiffusion within inhomogeneous media. For both models, in the 
presence of an external field, rare events induce a singular behavior in the originally Gaussian 
displacements distribution, giving rise to power-law tails. Remarkably, in the subdiffusive trap 
model, the combined effect of highly fluctuating waiting times and of a drift yields a non-Gaussian 
distribution characterized by long spatial tails and strong anomalous superdiffusion. 

PACS numbers: 05.40.fb 02.50.Ey 05.60.-k 



Large fluctuations and rare events play an important 
and often crucial role in many physical processes. Their 
contribution strongly influences physical obscrvables in 
many contexts, as for instance in laser cooling li- 
posome diffusion in ncmatic solutions 0, heteropoly- 
mers Q , delayed chaotic systems Q and non-equilibrium 
relaxation 

In standard diffusion, governed by the Fick's law, a 
Gaussian probability density function (PDF), P(x, t), de- 
scribes the position x of particles at time t, and the 
fluctuation-dissipation relations provide a powerful tool 
to investigate fluctuations by means of response measure- 
ments ■ In such a case, also the out-of-equilibrium re- 
sponse is well under control, as it is clear from the form 
of the PDF, which in the presence of an external field e 
reads: 
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where D is the diffusion coefficient. The field induces a 
finite drift, (x(t)) e /t = v e ^ 0, and a typical displacement 
xx{t), defined as the maximum of the PDF, which grows 
as Xt(1-) = v e t. 

Transport properties are different in the presence of 
large space inhomogeneities and/or several time scales, 
that can produce a deviation from the Gaussian diffu- 
sive behavior Large fluctuations in time can be ob- 
served in diffusion in strongly disordered systems, where 
particles can get trapped for a long time in a given po- 
sition. For instance in this paper we will consider the 
trap model which describes a particle diffusing with 
a power law distribution of waiting times, 



Pa(r) = T 
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tions characterize instead Levy-like motion in hetero- 
geneous and porous materials or in turbulent 
flow (l4^ , where transport is realized through increments 
of size I extracted from a power-law distribution p a (l) ~ 
H 1+Q > [Il-[l7|- 

0, if the tails 



and which is commonly studied as a prototype of slow 
dynamics in glassy systems 0. Large space fluctua- 



In the absence of an external field, e 
of the steps and waiting times distributions decay slowly 
enough, both the trap model and the Levy walk lead to 
anomalous transport, (x 2 (t)) ~ t 2 ^ z : subdiffusive, z > 2, 
for the trap model and superdiffusivc, z < 2, for the Levy 
walk. In these conditions, motion is largely influenced by 
rare events coming from the algebraic tails of the Pq{x, t), 
and, in the case where only one exponent is relevant, one 
can write the generalized scaling form: 

P (M)~^ 1/2 ^A 1/Z )- (3) 

Here F(y) is the scaling function and £(t) ~ t x / z is de- 
fined as the scaling length of the process. Assuming also a 
fast asymptotic decay for the tails of F(y), the anomalous 
behavior of the mean square displacement is straight- 
forward. However, we stress that, in general, F(y) can 
present a slow power-law decay, and the system can fea- 
ture a "strong" anomalous behavior, induced by the pres- 
ence of a physical cut-off in the distribution, which gives 
rise to a second scaling length in the system flil ]. 

When the fluctuations of the distributions are not too 
large and the tails decay sufficiently rapidly, Pq{x,£) is 
Gaussian and standard diffusion is recovered. Moreover, 
the standard equilibrium fluctuation-dissipation theo- 
rem, established in general also for anomalous dynam- 
ics holds with (x(i)) e growing linearly in time. 
One could then expect to observe a regular system in all 
its aspects and a modified PDF analogous to Q]) when 
an external field is applied. 

In this paper, we will show that the presence of a drift 
can induce an anomalous non- Gaussian dynamics even 
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for the regimes where equilibrium measurements show a 
standard diffusive behavior. In particular, by solving the 
master equations of the two processes mentioned above, 
we will determine the scaling forms of the PDF in the 
presence of an external field. As expected, we evidence 
that, in the region of anomalous transport, the distribu- 
tions are very sensitive to the presence of a drift. More 
surprisingly, in the regimes where the form of the distri- 
bution is simply Gaussian at equilibrium, we find that 
the presence of a field may trigger large fluctuations in- 
ducing a non-Gaussian shape of P t {x,t). This happens 
for I < a < 2 in the trap model and for 2 < a < 4 in 
the Levy walk: in these cases an arbitrary small external 
field e can induce a transition from standard to strong 
anomalous supcrdiffusion, namely to a regime in which 
the PDF and its higher order momenta are characterized 
by different scaling lengths. 

In the trap model, a particle moves with probability 
1/2 form x to x ± <5o, where, for simplicity, we take <5o 
constant. However, the main results do not change if 
the displacement follows a symmetric distribution with 
finite variance. Between two successive steps, the particle 
waits for a time r extracted from Eq. ([2]). In the Levy 
walk model again there are time intervals of duration 
r extracted from the PDF in Eq. ([2]), but in this case 
the particle, during each time lag, moves at a constant 
velocity vq, extracted from a symmetric distribution with 
finite variance, and hence performs displacements tvq. 

In the trap model, the external field is implemented by 
unbalancing the jump probability, i.e. setting to (l + e)/2 
and (I — e)/2 the probability of jumpin g to the right or 
to the left. For the Levy walk, in [22j, |23[ it has been 
shown that the most natural way to drive the system out 
of equilibrium is to apply an external field accelerating 
the particle during the walk, so that the distance traveled 
after a scattering event is 6 = vqt + er 2 . 

Inspired by the Gaussian case, the simplest generaliza- 
tion of Eq. (HJ) one may assume for e ^ reads as: 



P e {x,t) ~ r 1/z F[(x - v e t)/t 1/z }. 



(4) 



When e = 0, the left-right symmetry implies that v c = 
and F(y) is an even function. In standard diffusive pro- 
cesses, F{y) is Gaussian, z = 2, and the drift v e vanishes 
linearly with e. However, both in the trap model and in 
the Levy walk, Eq. ((4]) is made slightly more complex by 
the occurrence of rare events, which change the scaling 
structure introducing another characteristic length. Ac- 
cording to our results, the correct form for the PDF in 
the Levy walk is: 



P e (x, t) ~ r 1/z F[(x - v e t)/t 1/z ] Q (et 2 - x) 



(5) 



where O(x) is the Heaviside step function, which ac- 
counts for the maximum displacement possible up to time 
t, namely one single ballistic flight. In the case of the trap 
model the factor Q(et 2 — x) in Eq. ([5]) has to be replaced 
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FIG. 1: (color online) Collapse of the PDFs of the trap model, 
for q = 1.5 and field e = 0.15 according to Eq. (JJ) with 
velocity v t — 0.1. Symbols correspond to different times: 
t = 10 3 (•), t = 2 • 10 3 (■), t = 5 • 10 3 (♦), and t = 10 4 (A). 

Central inset: the same data of the main figure in log-log 
scale, as function of \x — v e t\. Notice the power law behavior 



P e (x, i) ~ x 



in agreement with the arguments reported 



in the text. Top left inset: collapse of the PDF for a = 1.5 
with zero drift, e = 0. Notice the simple Gaussian behavior, 
l{t) ~ t 1/2 , in agreement with Eq. ©. 



by Q(x), which cuts off the power law tail due to particles 
with large persistence time at the origin. Although these 
cuts do not affect the collapse of the displacement PDFs, 
as we will show, they are important in order to obtain 
the correct scaling of the moments (x n (t)) e , where (. . . ) e 
denotes the average over P e (x,t). 

The general scaling form for P e (x,t) in the trap model 
and in the Levy walk can be obtained from the master 
equation of the processes: this relates, at two subsequent 
scattering events, P e (x,t) with the function Q e {x : t), i.e. 
the probability of being scattered in x at time t [24| . In 
particular, the master equation for the trap model with 
vq extracted from a bimodal distribution, reads: 



<(x,t) 



ri + e 



1 - e 



Q e (x-S ,t-t') 



Q e (x + S 0> t-t') p a (t')dt 



+ 6(t)S(x) 



(6) 



with 



P e (x,t) 



(X,t-t') / Pa (T)dTdt', (7) 



where S(x) is the Dirac delta. The above equation can 
be solved asymptotically by a Fourier transform with re- 
spect to space and time, P(x,t) — > P(k,cj), in the equi- 
librium (e = 0) and out-of-equilibrium (e ^ 0) regimes 



Without drift, e = 0, the results are well known [24 1. 
In particular, in the asymptotic regime of large space 
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and time (k and u> close to zero), for the trap model one 
obtains: 



P ° {k,U) \(t)/(6 k 2 /2 + i(t)u>) if 



/(d k 2 /2 + C 2 Lu a ) if 0<a<l 



a > 1, 



(8) 

where, here and in the following, Cj and (t) are constants, 
determined by the distribution p a (r), and in particular 
(t p ) = jt p p a (t)dt. For a > 1 Eq. © corresponds to 
the Gaussian propagator of standard diffusion, while for 
< a < 1 the Gaussianity is broken and the dynamics is 
characterized by a subdiffusivc scaling length £{t) ~ t a l 2 . 

Let us now evaluate the propagators of the trap model 
in the presence of a drift, i.e. with e ^ 0. For large space 
and time we obtain: 

( dio^/iieSok + C 2 uJ a ) if < a < 1 
P e (fc,w) = < (i)/[i(e<5 fc + (*M + C 3 w a ] if Ka<2 
[(t)/[5 fc 2 /2 + i(e5 fc + if a>2. 

(9) 

Interestingly, all the regimes follow the general scaling 
pictures of equation (jU , however several different behav- 
iors are realized in the dynamics. For a > 2, the stan- 
dard response of a diffusive system is reproduced. In par- 
ticular, F(y) is a Gaussian moving at constant velocity 
v e = eSq/ (t). For < a < 1, where Gaussianity is broken 
already at equilibrium, F(y) is clearly non-Gaussian and 
the characteristic length grows as £(t) ~ t a . We notice 
that v e = 0, and hence a rigid motion of the probability 
with constant velocity is not observed. For 1/2 < a < 1, 
one finds a superdiffusive behaviour induced by the field, 
as discussed in In this case, due to the fast asymp- 
totic decay of F(y), the behaviour of momenta is ruled 
by £(t), namely (x n (t)) e ~ £ n (t), which is an example of 
weak anomalous diffusion (l8j . Comparing equations ([8]) 
and ©, for < a < 1 a singular behavior at e = is 
evident, i.e. F(y) changes its shape, as soon as £ / 0. 
In particular, the scaling function becomes asymmetric 
with respect to the origin y = 0, giving rise to an overall 
motion in the direction of the applied external field. 

The most intriguing case is for 1 < a < 2, where the 
presence of the field turns the diffusion from standard 
into strong anomalous: now the scaling function moves 
at constant velocity v e = eSo/{t), however the shape of 
F(y) is not Gaussian, and it develops a power-law tail 
on the left. This tail is responsible for the breaking 
of the relation (x n (t)) e ~ £ n (t) between momenta and 
scaling length of the PDF. In Fig. Q] we show the col- 
lapse of the PDFs at different times, according to the 
scaling for the case a = 1.5 with and without ex- 
ternal field (see inset). Care must be taken in doing 
the inverse Fourier transform of P e (k,ui). As suggested 
by numerical data and physical intuition, in real space 
F{y) must be cut off in y = and this provides a for- 
bidden zone also in the plane (fc,w) in calculating F(y) 
from P e (k,uj). Equation © implies then that the scal- 



ing length of the moving scaling function grows in a su- 
perdiffusive way, and in particular z = a < 2. The su- 
perdiffusive spreading of the PDF that one can see in 
Fig. [1] is triggered by the rare events corresponding to 
particles that at the observational time t are still at the 
origin. One can easily understand that the persistence 
of particles in the same position is directly related to 
the waiting time PDF ([2]). From this observation fol- 
lows that P e (x 7 t) ~ t\x — v e t\~( 1+a }0(x) (see central in- 
set of Fig. [T]). The cut 0(x) in the power law tail is 
necessary because at large times, even if particles have 
a non-negligible probability of staying close to the ori- 
gin, negative values of x are almost forbidden, due to 
the positive field. By considering the joint effect of the 
power-law behavior of the left tail and of the cut-off in 
zero, one can obtain the spreading around the peak of 



the PDF: ([5x(t)] 



(x 2 (t)) e - (x(t))l 



t 
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The 



strong anomalous behaviour is then clear because one 
has ([Sx(t)] 2 ) e ^£ 2 {t) ~t 2 ' a . 

Let us now consider the Levy walk. In this case the 
master equation is: 



Q e {x,t) = 



Q e (x-v t' -et' 2 ,t-t') 



+ -Q t {x + v t' -et' 2 ,t-t') p a {t')dt' 
+ 6{t)6{x) 



(10) 



with 

P e (x,t) 



rl 



Q e (x - v t' - et' 2 ,t - t 1 ) 



+ -Q t {x + v Q t' -et' 2 ,t-t') / Pa (T)dTdt' 



(11) 



At equilibrium, the Fourier transform Po(k,uj) reads: 



( ( 



C^-^iCak ' + C 6 oj a ) if < a < 1 
Po(k,u) = ( (t)/{C 7 k a +i(t)ui) if Ka<2 
[(t)/(v 2 (t 2 )k 2 /2 + i(t)io) if a>2. 

(12) 

In this case, standard Gaussian diffusion is recovered for 
a > 2, while for 1 < a < 2 and < a < 1 the non- 
Gaussian propagator corresponds to a superdiffusive and 
a ballistic motion, i.e. the scaling length grows as t 1 /" 
and as t, respectively. 

In the presence of a drift, from Eqs. (|T0|) and (|TT|) . in 
the asymptotic regimes where v € t <$^. t 2 / a we obtain for 
the Fourier transform P e (k,tj): 



P e (k,Gj) 



'C 8 u; Q ~7(C 9 fc Q / 2 + C 10 w Q ) if < a< 1 
(t) /(C 9 k a ' 2 + i{t)oj) if 1< a < 2 
(t) j\C % k a l 2 + i(e(t 2 )k + (t)u)] if 2 < a < 4 



,(t}/[(t 2 )v 2 k 2 /2 + l (e(t 2 )k 



(t)ui)] if a > 4 
(13) 
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From the above expressions one obtains the PDF in space 
and time by inverse Fourier transform, taking care also 
in this case of a forbidden region in the plane (k,ui). 
Indeed, at each observational time t, the largest dis- 
placement possible is x = vt + er 2 ~ er 2 , correspond- 
ing to the rare event of a particle never colliding up 
to r = t, and this bound necessarily requires for the 



G(et 2 



in the conjugated space must be taken with the con- 
straint |fc| > u 2 /e. With this caveat we find that for 
a > 4 Eq. (|13|) describes a Gaussian distribution, mov- 
ing at constant velocity v e = e(t 2 )/(t). For a < 2, the 
rigid translation of the scaling function is subdominant 
(v t t <C t 2 / a ) and the global effect of the field is a strong 
asymmetry of F(y) (see Fig. [2]). In this case the scaling 
length grows in super-ballistic way, i.e. z = 1/2 (lin- 
early accelerated motion) for < a < 1 and z = a/2 for 
1 < a < 2. Finally, for 2 < a < 4 the drift with con- 
stant velocity v e becomes relevant, but the scaling func- 
tion is still not Gaussian due to an emerging power-law 
tail which produces a strong anomalous superdiffusive 
behavior. The remarkable feature we discovered for the 
Levy walk with 2 < a < 4 is that, similarly to what hap- 
pens for the trap model with 1 < a < 2, the field of in- 
tensity e produces a power-law tail otherwise not present 
in the zero- field Gaussian distribution (see Fig.|3]). More- 
over the scaling lenght £(t) which yields the collapse of 
the PDFs does not govern the behaviour of momenta, 
(x n (t)) e l n {t). One can obtain that the probabil- 
ity of the largest displacement x at time t goes like 
~ x _ ( 1+a / 2 ) and such a power-law behavior of the tail can 
be clearly seen in Figs. [2] and [3] Indeed this power-law 
behavior, together with the cut-off Q(et 2 — x), provides 
for the spreading around the peak of the distribution: 

([Sx(t)] 2 ) e = (x 2 (t)) e - (x(t)) 2 ~ t 5 - * e(t) 2 ~ t^, 

showing what claimed above: a single scaling length can- 
not capture the behavior of all momenta. 

In this letter we have presented results for the trap 
model and the Levy walk in the presence of currents. 
At variance with standard Gaussian systems, in such 
models the external field, even arbitrarily small, signifi- 
cantly modifies the scaling properties of the PDF of dis- 
placements, introducing a new lengthscale related to rare 
events. Beyond the principal scaling length of the dis- 
tribution £(t) and the one related to the rigid shift (if 
present), one must also consider the typical length in- 
troduced by the cut-off of the power-law tail, necessary 
for the calculation of higher order momenta. That is 
how rare events induce the observed strong anomalous 
behaviour. 

In conclusion we have shown that the modification of 
transport properties in the presence of an external field 
represents a valuable probe to unveil the underlying dy- 
namical structure of the system. In particular, field- 
induced anomalous behaviors highlight the importance 
of rare and large fluctuations in regimes where the diffu- 
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FIG. 2: (color online) Collapse of PDFs with the scaling 
t 2 / a P t (x/t 2 / a ), for the Levy walk with exponent a = 1.5 and 
field e = 0.25. Inset: P e (x,t) at different times. Notice that 
the approximation (x — v e t)/t 2 ^ a ~ x/t 2 ^ a works very well 
asymptotically. 
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FIG. 3: (color online) Collapse of the PDFs according to the 
scaling t 2/a P t [(x — v t t)/t 2/a ], for the Levy walk with expo- 
nent a = 2.5 and field e = 0.25, with v t = 0.38. Symbols 
correspond to different times: t = 1.6 ■ 10 4 (•), t = 3.2 ■ 10 4 
(■), t = 6.4 ■ 10 4 (♦), and t = 1.28 ■ 10 5 (A). Inset: the 
same data of the main figure are reported in log-log scale to 
highlight the power-law tail of the distribution. 



sional properties are apparently standard. 
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